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This study proposes a thermoelectric cooling device based on an atomic-sized junction. Using first- 
principles approaches, the working conditions and coefficient of performance (COP) of an atomic 
junction as an electronic refrigerator are investigated. Quantum nature due to the size minimization 
offered by atomic-level control of properties facilitates electron cooling beyond the expectation of the 
conventional thermoelectric device theory. Our studies reveal that the 4-Al junction with ZT=0.02 

■ could have refrigeration efficiency comparable to a bulk thermoelectric refrigerator with a ZT of an 
(N . order of magnitude larger. 
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I. INTRODUCTION 

' The transport properties of atomic-scale junctions are curently eliciting a tremendous wave of multidisciplinary 
f~| ', scientific interest [l|, |2| . This interest is motivated by the aspiration to develop new forms of electronic devices based 
^ • on subminiature structures and to understand the fundamental properties of non-equilibrium transport in atomic- 
Qj [ scale junctions. A growing number of studies are currently available to diversify the scopes of transport properties in 
molecular /atomic junctions, including current-voltage characteristics fsUvIl, inelastic electron tunneling spectroscopy 
. (lETS) [8-lii] , shot noise [l7l-l20|. counting statistics local heating_|22,,^,23j. and gate-controlled effects [24-28]. 
, Substantial progress in experiment and theory has also been achieved [29l - l3l| . 

■ Recently, considerable attention is being paid to thermopower of nanojunctions. Atomic scale thermoelectric 
I , junctions are gaining increased attention due to the recent measurements of Seebeck coefhcient, defined as = dV/dT, 

' in metallic atomic/molecular junctions [32| - [35j . Measurements of Seebeck coefficient provide a useful experimental 
Ch approach to exploring the electronic structure of the molecule bridging the electrodes |3a|. These experiments inspire 
O rapid development in the theory of thermoelectricity at the atomic and molecular scale |37l-[49| . 

I I Exploring the possibility of subminiature energy-conversion devices is scientifically and technologically important. 

Thermoelectric effects hybridize the electron and energy transport, which complicates the fundamental understanding 
for quantum transport of electron and energy under non-equilibrium conditions. Thermoelectric nanojunctions open 
^ new challenges for developing novel energy-conversion nano-devices, such as thermoelectric power generators and 
refrigerators at the atomic/molecular level. Methodologically, the scope of the research needs to extend through the 
Q.^ [ utilization of unprecedented experiment and theory. Thus far, previous reports, focus only on the Seebeck coefficient 
1 (S) and the thermoelectric figure of merit (ZT) of nanojunctions. Little attempt is made to explore new forms of 
thermoelectric devices based on atomic-scale junctions i50, SlJ- This project initiates research on the theory of working 
conditions and coefficient of performance (COP)in thermoelectric nano-refrigerator. It offers a key to understanding 
the nanoscale energy-conversion mechanism and may open a new path in exploring novel thermoelectric nano-devices. 
, The Seebeck coefficient is related closely to the magnitude and slope of transmission function near the chemical 
' potentials. To have an efficient energy-conversion device, the thermoelectric junctions with large magnitudes of 
Seebeck coefficients and small magnitudes of thermal conductances are of key importance. Such systems are often 
characterized by sharp density of states (DOSs) near the chemical potentials and poor mechanical coupling between 
^ I the nano-object and electrodes. 

Compared with atomic-scale thermoelectric power generators [52j . thermoelectric refrigerators require more strin- 
gent working conditions. In the bulk and mesoscopic system, the efficiency of a thermoelectric (TE) refrigerator is 
usually suppressed by the large work function and the overwhelmingly irreversible Joule heating. For example, the 
operation of a thermoelectric cooling device, such as a vacuum diode, is limited to very high temperatures due to its 
large potential barrier [53| . Recently, researchers have shown the prospect of operating thermoelectric refrigerators 
at room temperature using semiconductor heterojunctions to reduce the work function [s^ - fSTj . An alternative new 
approach to low temperature-operated thermoelectric refrigerators may be atomic-sized junctions, the extreme limit 
of device miniaturization. 
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This study addresses the question of thermoelectric coohng based on an atomic junction operated at the apphed 
bias Vb < Vonset, where Vonset is the onset bias of local heating. To enrich the understanding of thermoelectric 
cooling mechanism at the atomic level, we develop an analytical theory which offers more elaborate explanation of the 
working conditions, threshold bias, and behaviors of COP in thermoelectric nano-refrigerators. As an example, we 
investigate the thermoelectric cooling effects in the 4-Al junction as an electronic refrigerator with and without the 
thermal current carried by phonons using first-principles approaches. The atomic thermoelectric junction is marked 
by the reduced tunneling potential due to bridging atoms, the enhanced Seebeck coefficient due to a Px-Py sigma 
channel, the suppressed local heating due to discrete vibrational modes, and the reduced photon radiation due to 
small contact surface. Due to size minimization, these quantum natures remarkably facilitate efficient thermoelectric 
cooling effect beyond the expectations of the conventional thermoelectric device theory. Our calculations indicate 
that the novel thermoelectric nano-refrigerator is promising because it allows electronic nano-refrigerators to perform 
in a wide range of temperatures, even at temperatures under 100 K, which is a great improvement in comparison with 
the vacuum diode [ssj . 

The cooling effect in thermoelectric nanojunctions could be destroyed by the adverse phonon's thermal current, 
which is likely to depend strongly on the detailed atomic structure of the full system especially in the contact region. 
To minimize the adverse effects imposed by the phonon's thermal current, we suggest creating a weak mechanical link 
between the nano-structured object and the electrodes while still allowing electrons to tunnel. Our theory reveals 
that the choice of a thermoelectric junction with a large Seebeck coefficient is of crucial importance to the design of 
nano-refrigerator. Widely diversified atomic-sized junctions may be achieved by manipulating the species of nano- 
structured objects and contact regions. Such manipulations may lead to a significant change in DOSs, consequently 
varying the Seebeck coefficients of nano-j unctions. Atomic-level control of the contact region plays a primary role in 
yielding enhanced Seebeck coefficient and creating poor thermal contact. A full exploration of all possibilities in such 
an unknown system may open new opportunities and challenges in developing new forms of thermoelectric cooling 
devices at the atomic level. 

The flow of the discussion in this paper is as follows. Firstly, we describe the details of density functional theory, 
theory of thermoelectricity, and theory of thermoelectric nano-refrigerator in the absence and presence of phonon's 
thermal current in section II. Theoretical Methods. In section III. Results and discussion, we discuss the 
thermoelectric properties of 4-Al atomic junctions, and subsequently discuss the onset bias, working conditions, and 
COP of thermoelectric nano-refrigerators. We then summarize our findings in section IV. Conclusions. 

II. THEORETICAL METHODS 

The theories presented in the following subsections are general to any atomic/molecular junction that utilizes 
tunneling as the major transport mechanism, provided that the applied bias is much smaller than the onset bias 
of local heating. In the subsection A. Density Functional Theory, we present an introduction of the density- 
fimctional theory (DFT). In the subsection B. Theory of Thermoelectricity, we present the theory to calculate 
the Seebeck coefficient, electric conductance, electron's thermal conductance, phonon's thermal conductance, and 
thermoelectric figure of merit {ZT). In the subsection C. Theory of Thermoelectric Refrigerator, we present 
the theory for the rate of thermal energy removal from the cold-temperature reservoir, the working condition for 
thermoelectric refrigeration, and the coefficient of performance for the atomic-scale thermoelectric refrigerator with 
and without the presence of thermal current of the phonon. 

A. Density Functional Theory 

The electric current, electron's thermal current, and Seebeck coefficient are calculated in a truly atomic-scale 
junction in terms of the effective single-particle wave-functions obtained self-consistently within the DFT. In this 
subsection, we present a brief introduction of how to calculate the electric current and the electron's thermal current 
carried by the electron transport in the DFT framework. We picture a nanoscale junction as formed by two semi- 
infinite electrodes held apart at a fixed distance with a nano-structured object bridging the gap between them. The 
full Hamiltonian of the system is H = Hq + V, wherein Hq is the Hamiltonian due to the bare electrodes and V is 
the scattering potential of the nano-structured object. 

First, we calculate the wave functions of the bare electrodes with an applied bias Vb = {fJ^R — /iL)/e, where fJ.L{R) is 
the chemical potential deep in the left (right) electrode. The unperturbed wave functions of the bare electrodes have 
the form, vl/^^(^)(r) = e«K„.R 

■""bk^ii (^)' where u^'--^j{z) describes the electrons incident from the left (right) electrode 
before the inclusion of the nano-structured object. The equation u^^^\z) is calculated by solving the Shrodinger 
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FIG. 1: (Color online) (a) Scheme of the 4 — Al atomic junction with different temperatures Tl{r) and chemical potentials 
t^L(R)- The Al-jellium distance is about 2.5 a.u. and the Al — Al bond distance is about 6.3 a.u.; (b) Density of state and 
projected density of states of the i — Al junction; (c) Seebeck coefficient vs. temperature of the 4 — Al junction; and (d) ZT 
vs. T with the phonon's thermal current for various values of A. 



equation and Poisson equations iteratively until self-consistency is obtained. Note that u^-^^^ (z) satisfies the following 
boundary condition: 

where Ky is the electron momentum in the plane parallel to the electrode surfaces, and z is the coordinate parallel 
to the direction of the current. The condition of energy conservation gives ^/c^ = E — ^ l-'^lll^ ~ ''^e//(oo) and 
— E — ^ l-'^lll^ ~ ''^e//(~oo), where Vf,ff{z) is the effective potential comprising the electrostatic and exchange- 
correlation potentials. 

The nano-structured object is considered in the scattering approaches. The scattering wave functions of the entire 
system are calculated by solving the Lippmann-Schwinger equation iteratively until self-consistency is obtained, 

*^K,!W = ^'^^fir) + J d'r, J d'r,G%{r,r,)Vir,,r,)^'^<-^lir,), (2) 

where ^EK^li^) stands for the effective single-particle wave functions of the entire system, which also represents the 
electrons with energy E incident from the left (right) electrode. The potential y(ri,r2) that the electrons experience 
when they scatter through the nanojunction is 

V{r,,r2) = Vps{ri,r2) + Uv^c [n (ri)] - V^c K (ri)]) + J dv3-^^^^\ S{ri - rs), (3) 
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where V^s(ri, is the electron-ion interaction potential represented with pseudopotential; Vxc [n (r)] is the exchange- 
correlation potential calculated at the level of the local-density approximation; ng (r) is the electron density for the 
pair of biased bare electrodes; n (r) is the electron density for the total system; and Sn (r) is their difference. The 
quantity is the Green's function for the bare electrodes. The wavefunctions that achieve self-consistency in the 
DFT framework in plane wave basis are applied to calculate the electric current, and the electron's thermal current. 

These right- and left-moving wave functions weighed with the Fermi-Dirac distribution function according to their 
energies and temperatures are applied to calculate the electric current as 



eh 



dE / dY±_ / dK 



LL 



(4) 



where Iee''^^1{y) 



,R{L) 



(r) 



,T,i?,(L) / X 



^"^/'^^ (r) and drj_ represents an element of the 



electrode surface. Here, we assume that the left and right electrodes are independent electron reservoirs, 
with the population of the electron described by the Fermi-Dirac distribution function, /^(/^^(i?)! 2^L(fl)) — 
1 / {e:xp[(^E — ^E(R)) / (f^BTi^iR))] + 1}, where Hl{r) and Ti^{r) are the chemical potential and the ternperature in 
the left (right) electrode, respectively. More detailed descriptions of theory can be found in Refs. [3. l65j|. 
The above expression can be cast in a Landauer-Biittiker formalism: 



2e f 

I{fiL,TL;HR,TR) = / dET{E)[fE{fiR,TR) ~ fE{fJ.L,TL)], 



(5) 



where t{E) = t^{E) = t^{E) is a direct consequence of the time-reversal symmetry, and t^'^'"\E) is the transmission 
function of the electrons with energy E incident from the right (left) electrode. 



R(L) 



{E) 



J dr^J dK,|/-f-^(r). 



(6) 



There is an analog between the electric current and electron's thermal current. The flow of elections can also transport 
energy. The electron's thermal current, defined as the rate at which thermal energy flows from the right (into the 
left) electrode, is 



Jel'^\tJ'L,TL]lMR,TR) = ^ J dE{E - ^Lt^ii))T{E)[fE{HR,TR) - fE{tJ-L,TL)]. 



(7) 



B. Theory of Thermoelectricity 

In this subsection, we present the theory to calculate the zero-bias electric conductance, Seebeck coefHcient, and 
electron's thermal conductance for a thermoelectric nanojunction. These quantities have been calculated in terms 
of the effective single-particle wave-functions obtained self-consistently within the static density-functional theory in 
truly atomic-scale junction. 

We assume that the left and right electrodes serve as independent temperature reservoirs. The population of 
electrons in the left (right) electrode is described by the Fermi-Dirac distribution function, fE{fJ'L(R},TL{R)) — 
l/{exp[(i? — fJ-L(R.)) / {kBTi^{R))] + 1}, where the chemical potentials are fiR = fiL = fJ- (i.e., at zero external bias). Let 
us now consider an extra infinitesimal current induced by an additional infinitesimal temperature (dT) and voltage 
(dV) across the junctions in a open circuit. The current induced by dT and dV are {dI)T = I{^J', T; fi,T + dT) and 
{dl)v = I{fJ,,T;fi + dV,T), respectively, where I{n,T; fi,T + dT) and I{id,T; i^,T + dT) are given by Eq. Suppose 
that the current cannot actually flow in an open circuit, thus, {dI)T counterbalances {dl)v- In other words, the extra 
net current is zero, 

dl = /(^, T;fi + edV, T + dT) ^ {dI)T + {dl)v ^ 0. (8) 

The Fermi-Dirac distribution function in Eq. ([5]) can be expanded up to the first order in dT and dV , and we obtain 
the Seebeck coefficient (defined by = dV/ dT) , 



where 



K„i^^,T) = - / dEr{E) {E - ^)" (10) 
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The Seebeck coefScient in the low-temperature regime can be obtained by expanding Kn{n,T) to the lowest 
order in temperatures through the Sommerfeld expansion, that is, Kq « ''"(/i), Ki « [7r^fc^T'(/i)/3]T^, and 
K2 ~ [tt^ kgT{iJ,) /S\T^ . The Seebeck coefficient up to the lowest order in temperature is 

S « aT, (11) 

where a = —n^kg / (3er(^)). The Seebeck coefBcient is positive (negative) when the slope of transmission 
function is negative (positive), which is closely related to the transmission function near the chemical potentials. We 
note that the Seebeck coefficient can be optimized via the gate voltage which shifts the DOSs. 

The electron's thermal current is the energy current carried by electrons traveling between electrodes driven by dT 
and dV . Analogous to the extra current given by Eq. ([S]), the extra electron's thermal current is 

dJel = idJel)T + {dJ,i)v. (12) 

where (dJei)T = Jel{^^^T] fj,,T + dT) and {dJei)v — Jeiil^TT; fi + edV,T) are the fractions of electron's thermal current 
driven by dT and dV, respectively. Note that dV is generated by the Seebeck effect according to the temperature 
difference dT. Both Jei {n,T; ii,T + dT) and Jei {n, T]^ + edV, T) can be calculated using Eq. (O . 

Given that we define the electron's thermal conductance as fcg/ — dJei/dT, the electron's thermal conductance k^i 
can decomposed into two components: 

Keiiii,T) = K^i{fi,T) + h;^i{^i,T), (13) 

where — {dJei)T / dT and = {dJei)v/dT. We note that kJj and are the portions of the electron's thermal 
conductance driven by dT and dV, respectively. Analogous to Eq. ([9]), k^i and can be can be expressed by Eq. (|10p : 

J, 2 K2{ti,T) 

l^elif^^T) = - , (14) 

and 

K^li^l,T) = ^K,i^,,T)S{f,,T). (15) 

One should note that = if the Seebeck coefficient of the system is zero because dV is zero. 

In the low-temperature regime, kJ^ and can be expanded to the lowest order in temperatures using the Som- 
merfeld expansion: 

« /3vT ^and w I3tT, (16) 

where Pv = — 27r^fc^[T'(/u)]^/(9/iT(^)) and Pt = 27r^fc^T(^)/(3/i). In the above expansions, we also applied the 
following approximations: Ki{^i,T) « [Tr^k%T'{fi)/3]T'^, K2{n,T) « [7r2fc|r(^)/3]T2, and Eq. (HH). In the low- 
temperature regime, dominates the electron's thermal current. Thus k^i w is linear in T, that is, 

fieiit^,T) ^ I3tT. (17) 

At zero bias, the electric conductance can be expressed as 

In the low-temperature regime, the zero-bias conductance is usually insensitive to temperatures if tunneling is the 
major transport mechanism. 

Thus far, the physical quantities (S*, a, and k^i) previously discussed are related to the propagation of electrons. 
The heat current carried by phonon may occur in a real system. The phonon's thermal current, which is driven 
by the temperature difference AT, flows from the hot reservoir into the cold reservoir. To determine the impact of 
the phonon's thermal current on refrigeration, the weak link model is chosen to describe it. The weak link model 
assumes that the nanojunction is a weak elastic link with a given stiffness K that can be evaluated from total energy 
calculations or from experimental measurement |6(i]. Two metal electrodes are regarded as the macroscopic bodies 
under their thermodynamic equilibrium, and are taken as ideal thermal conductors. To the leading order in the 
strength of the weak link, the phonon's thermal current {Jph) via elastic phonon scattering is, 

Jph = J dEENL{E)Nn{E)[nL{E) - nR{E)], (19) 
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where K = AY /I is the stiffness of the weak elastic hnk; Y is the Young's modulus of the junction; A is the cross- 
section; I is the length of the nano-structured object as the weak elastic link bridging the electrodes; Ni^i^ji^ (E) is 
the spectral density of phonon states at the left (right) electrode surface which is measurable by experiments; and 
n^(fl) = l/(e^/^^"^'^(«) — 1) is the Bose-Einstein distribution function. In the long wave length limit, the spectral 
density of surface phonon is given by Ni^(^i^^{E) « CE. The phonon's thermal conductance defined as Kp^ — AJph/AT 
as AT 0, is 

AT AT 

Kph = [Jph {Tl + —,Tr-—)~ Jph {Tl ,Tr)]/AT. (20) 

Expanding the Bose-Einstein distribution function in Jph{TL + T^ — ^) to the first order of AT, we obtain the 
phonon's thermal conductance: 

When TrkTl^ T, then Eq. (HH) is reduced to, 

= dEE^^. (22) 

where n = l/[e'^/^''BT) _ 

The efficiency of thermoelectric nano-devices is conventionally described by the thermoelectric figure of merit ZT , 
which depends on the following physical factors: Seebeck coefficient (S"), electric conductance (cr), electron's thermal 
conductance (Kgj), and phonon's thermal conductance (k^/i). ZT is defined as 

ZT = — T, (23) 

where S, a, Kei, and Kph can be numerically calculated using Eqs. ([9]), (ITSl) . (fT3)) . and ([22]). respectively. When ZT 
tends to infinity, the thermoelectric efficiency of nanojunctions will reach Carnot efficiency. 



C. Theory of Thermoelectric Refrigerator 



In this subsection, we present the thermoelectric cooling device theory at atomic scale. The left- and right-electrodes 
are assumed to be independent electron reservoirs, respectively, with the electron's population described by the Fermi- 
Dirac distribution function, respectively. The left (right) electrode serves as the hot (cold) temperature reservoir with 
temperature Tr = Tq (Tl = Th = Tc + AT) and the phonon's population is described by Bose-Einstein distribution 
function. We consider the nanojunction connecting to an external battery with bias Vb = {^J■R — fJ.L)/e, which drives 
the electrons fiowing from the right- to left-electrodes. The thermal current carried by the electron traveling between 
two electrodes is given by Eq. It should be noted that 



jSifJ-L,TL;HR.,TR) + I{fiL,TL;fiR,TR)VB = J^i{fiL,TL;fiR,TR), 



(24) 



which is due to the energy conservation, where I{^l,Tl; ^iR, Tr) is the current given by Eq. ([5]). Equation (|24|) states 
that it consumes electric energy to take heat from the right (cold) into left (hot) reservoir. Thus, the thermoelectric 
junction can be regarded as an electronic cooling device when > 0, which states that the thermoelectric junction 
is capable of removing heat from the cold reservoir. 

A measure of a refrigerator's performance is the ratio of the rate of heat removed from the cold reservoir to the 
electric power done on the system. The ratio is called the coefficient of performance (COP): 

jR jR 

(25) 



where Vb is the bias applied across the nanojunction, / is the electric current, and is the rate of thermal energy 
removed from the cold reservoir. The thermoelectric junction as a nano-refrigerator is working when > (and 
thus ry > 0). 
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1. Properties of Thermoelectric Refrigerator in the Absence of the Phonon's Thermal Current 



In the following, we develop an analytical theory to gain insight into the fundamentals of the cooling effect in the 
thermoelectric nanojunction. We apply the Sommerfeld expansion to Jf^(AiL, 7l; /i^, Tr) and obtain, 

- I [r(/iL)(Mfi - Mi)'] - ^ [r'{f^L){m - ml)'] , (26) 

where we use the following relations: t{E) t{^l)+t' [hl)[E- ^l), t{iib:) ^ t{^l)+t' {^iL){^J-B- fJ-L), Tl = Tr + AT, 
Jl^^iE - fiR)dE = ~^{^iR - jl^^{E - ^lL){E - m)dE = -\{jiR - ^iL)^ and Mfl - Ml = eVs. When AT « T^, 

Equation can be expressed as a polynomial of Vb '■ 

J^i=-a + bVs - c^l - dVi , (27) 

where a = 2Tr'^ k%T{n)Tn AT / {3h), b = 2e7r2fc|T'(^)r|/(3/i), c = e^T(p)/h, and d = e3T'(/i)/(3/i). The above equation 
is convenient for analytical exploration of the properties of thermoelectric nano-refrigerators. 

A thermoelectric nano-refrigerator functions only when the maximum value of the electron's thermal current is 
positive, that is, {J^^)^^^^ > 0. In a small bias regime, the term Vg can be neglected in Eq. (P7)) . i.e., k, 
—a + bVB — cVg. In this case, the working condition of nano-refrigerator is given by {Jei)ina.^ ~ (— 4ac + b^) / (4c) > 0, 
which yields the criterion for the existence of electronic cooling 



where S = — ^^J^ \('^) '^R Seebeck coefficient of the nanoscale junction [isj . 

For a given Tr and a given temperature difference AT, J^{Vb) is a function of bias Vb- We observe that is 
negative at Vb = 0, and thus a lower limit of bias (denoted as V^^') is needed. The threshold bias is defined as the 
smallest positive solution of J^{Vb) = 0. For Vb > V^^, > and thermoelectric cooling effect exists. V^f has a 
very small value; hence, we keep the terms up to 0{Vb) in Eq. ([27|). i.e., J^{Vb) w — a + 6Vb. In this simple case, 
the threshold bias V^J^ w a/b is derived from J^i{V^l) = 0, from which we obtain, 

V,t = -^AT, (29) 

where we have assumed AT ^ Tr. Eq. p9| shows that the threshold bias V^^ is independent of Tr and increases as 
AT increases. 

We note that is a function of bias Vb, Tr, and AT. For a given bias Vb and a given temperature difference AT, 
is a function of Tr. We observe that there is a lower limit of temperatures when the thermoelectric refrigerator 
is working. The onset temperature for refrigeration effect, denoted as Top, is defined by J^(Top) — 0. The nano- 
refrigerator is not functioning {J^{Tr) < 0) when Tr < Top. Especially, J^{Tr) can be expressed as a polynomial of 
Tr, derived from Eq. If the Seebeck coefficient is sufficiently large and we neglect the terms higher than O(T^), 

the threshold operation temperature Top can be calculated analytically by solving the polynomial J^{Tr) = 0, which 
gives: 



Top = HVb + VWVbW+Wb] , (30) 



where a — and /3 ~ jr^k'^i'^lfj,) ■ Equation (|30|) shows that Top increases as AT and Vb increase, respectively. 

When AT = 0, a = and Eq. 1^ approaches to 



VbTr 
2S 



('^op)min ~ \/^7Er- (31) 



Equation (PT|) shows that Top approaches to the lower limit {Top)^-^^ as AT — > 0, where {Top)^-^^, increases as Vb ■ Tr 
increases. 
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We now turn to investigate the COP of the thermoelectric coohng device. When Vb is smah, jj^ — = IVb ~ crVg 
. In this case, the COP is in the fohowing form, 

- ^> (32) 

where JJ'(Vb) is a polynomial of Vb given by Eq. (P5|) . The applied biases considered in this study are small, hence 
we can consider J^{Vb) [given by Eq. p7)) ] up to 0{Vb)- Thus, Eq. (|32|) takes the form, 



"•■^--^^ ,33, 



The maximum value of COP (denoted as ?7l^ax) occurs at ?7'(V^max) = 0; which gives 

W 4^2fc|AT 2' ^ ^ 

where the maximum value of 77*^' occurs at bias Vb — V^axj where 

^Zx « 214t (35) 

where V^'^' is given by Eq. ([29)) . Concluding from Eqs. ((29| and p4l) . the refrigeration effect is triggered at around 
Vb = Vff and optimized around Vb = SVj.^'. 



2. Properties of Thermoelectric Refrigerator m the Presence of the Phonon's Thermal Current 

Phonon's thermal current is an adverse effect to thermoelectric refrigeration because it heats up cold electrodes. 
To realize the impact of this adverse effect to refrigeration, we consider the weak-link model suitable for describing 
the heat transport for two thermal reservoirs connected by a weak elastic link 60]. In the low-temperature regime, 
Eq. (j2ip can be expanded up to the lowest order in temperatures, 

= A(T| - Tt), (36) 

where A ~ 2tt^ K'^ C'^ k% / {15fi) , where K is the stiffness of the nano-structured object bridging the metal electrodes, 
and C is the slope of the surface phonon's dispersion function in long wave-length limit. The effect of the phonon's 
thermal current on refrigeration is described by a single parameter A, which is determined by K and C. 

Observed from the recent experiments [62 ] , the stiffness K of the linear Pt atomic chains varies in a wide range of 
several orders of magnitudes. The reason for this experimental finding could be that the atomic chain is particularly 
stiff along the chain direction in a perfectly linear atomic chain; otherwise, the atomic chain could be weak and 
vulnerable to bending. The magnitude of phonon's thermal current is sensible to the detailed atomic geometries in 
the contact region and could be uncertain. The simplified weak link model allows us to use a single parameter A to 
investigate the effect of phonon's thermal current on refrigeration. 

To determine the impact of the phonon's thermal current (Jj^), we repeat similar discussions in the previous 
subsection for the thermoelectric nano-refrigerators. Correspondingly, the COP of the thermoelectric nano- refrigerator 
becomes 

jR 

IVb 

where J^^^i^ — ^ei ^ ^ph combined thermal current. 

We assume that AT ^ Tr and limit ourself to small bias regime {Vb <?C 20 mV) to avoid substantial local heating. 
In this case, Eq. (|36l) is reduced to. 



J^,, « -4ArlAr. (38) 

Using Eqs. pT]) and ([38]). the combined thermal current can be expressed in terms of a polynomial of Vb similar to 
Eq. (EZD, 

= -aei+ph + bVB - cVl - dVl, (39) 
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where aei+ph = 2Tr'^k%T(^i)TRAT/{3h) - 4AT|AT; the coefficients b, c, and d remain the same as those in Eq. (j27|) . 

An analogy can be drawn here between Eqs. ([5^ and (P7|) . and Ukewise between Eqs. ([57]) and Consequently, 
the discussions in the previous section can be easily replicated here. In the presence of phonon's thermal current, the 
working condition of the nano-refrigerator becomes, 



•2^ 8r|\ / AT 



where A represents the strength of the phonon's thermal current, and a is the conductance of the nanojunction. As 
the phonon's thermal current vanishes (A = 0), Eq. (HUl) restores Eq. Equation (PU]) shows that the phonon's 

thermal current is an adverse effect to thermoelectric refrigeration. A sufficiently large phonon's thermal current has 
a large value of A, which is likely to break down the inequality described in Eq. (PO)) and ruin refrigeration capability. 
For a given Tr and a given temperature difference AT, Jq {Vb) is a function of bias Vb- We note that Jq is negative 

at Vb — 0; therefore, a minimum bias (denoted as V^^^'^'^'^) is needed to trigger possible thermoelectric cooling effect. 
The threshold bias V^l'^^^ is defined as the first zero of Ji^(V^B) = 0. Owing to small value of V^|^'^^'' and AT < Tr, 
the threshold bias becomes 

Vtf" = Vtt + A [4 i-S) Tl^T/a] , (41) 

where V^^ is given by Eq. ([29]) which is analogous to Eq. (|41|) . The threshold bias V^l^^^^ increases as the intensity 
of the phonon's thermal current (A) increases. 

Finally, we investigate the optimized value of the COP for a given Tr and a given temperature difference AT. The 
upper limit of COP considering the phonon's thermal current is given by 



el+ph _ 
/max 



where ?7max is the upper limit of COP given by Eq. (p4)) which is analogous to Eq. p2)) . 



(42) 



III. RESULTS AND DISCUSSION 



The aim of this research is to investigate the novel thermoelectric cooling effect based on a thermoelectric junction 
at the atomic level. First-principles calculations and analytical theories are proposed in order to apply the quantum 
natures offered by the nano-structures to the function of thermoelectric nano-refrigerators. The atomic-scale ther- 
moelectric junction is marked by the enhanced Seebeck coefficient due to reduced tunneling potential via bridging 
atoms, the suppressed local heating due to the nanoscale structure as well as small operation bias, and the suppressed 
photon radiation due to small size. Atomic-level control of the contact region could provide an alternative method 
to suppress the phonon's thermal current via weakened mechanical coupling between nano-structured object and 
electrodes. The above notable features suggest that the atomic-scale thermoelectric junctions as nano-refrigerators 
are worth examining. 



A. Thermoelectric Properties of the 4-AI Atomic Junction 

We propose a 4-Al junction as a thermoelectric cooling device, as depicted in Fig. [IJa). The aluminum junction 
is an ideal test bed for comparing quantum transport theories under non-equilibrium and experiments (ssl , |59| . The 
proposed 4-Al junction has a salient Seebeck effect due to the Px-Py orbital near the chemical potentials. The 
4-Al thermoelectric junction serves as an examples illustrating the advantage of nano-scale refrigerators, where the 
overwhelmingly Joule heat in the bulk system and photon radiation are strongly suppressed due to size reduction. The 
phonon's thermal current, which is relevant to the mechanical elasticity of nanojunction, could also be suppressed by 
controlling the atomic contact region. These properties facilitate possible thermoelectric cooling in the 4-Al junction 
beyond the expectation of conventional solid-state device theory. 

We begin our discussion by considering an ideal 4-Al atomic chain bridging two bulk Al metal electrodes that we 
model as ideal metals (jellium model, Vg « 2). The nano-structured object is considered in scattering approaches. The 
scattering wave functions of the whole system are calculated by solving the Lippmann-Schwinger equation iteratively 
until self-consistency is obtained. We assume that the left and right electrodes are independent electron and phonon 
reservoirs (Tc = Tr; Th = Tl = Tr -\- AT), respectively, with the electron and phonon population described 




FIG. 2: (Color online) (a) Critical operation temperature Top vs. AT and Vb in the absence of the phonon's thermal current; 
(b) Top vs. AT and Vb in the presence of the phonon's thermal current with A = 10~^^ Watt/K'*. 
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by the Fermi-Dirac and Bose-Einstein distribution function. Two electrodes are connected to a battery with an 
apphed bias Vb — (/ifl — fJ^L)/^, where hl and fin are the chemical potential deep in the left and right electrodes, 
respectively.A detailed account of the theory is given in subsection II. A Density Functional Theory. The choice of 
thermoelectric junction with a large Seebeck coefficient is of crucial importance for the design of thermoelectricnano- 
devices. Therefore, we investigate the thermoelectric properties of the 4-Al atomic junction as follows. 

The Seebeck coefficients are calculated from first-principles using the transmission function obtained from the DFT 
calculations, as described in Eq. The 4-Al atomic junction is marked by a sigma channel near the chemical 
potentials with the Px-Py orbital characters, as shown in Fig.[IJb). The Seebeck coefficients correlate highly with the 
magnitudes and slopes of DOSs near the chemical potentials, as described in Eq. pT|) . The sigma channel leads to a 
large value in the slope of the transmission function near the chemical potentials. This is the reason for the enhanced 
Seebeck coefficient in the 4-Al atomic junction, as shown in Fig.[TJc), which compares favorably with those of Pt, Pd, 
and Au atomic chains. 

The efficiency of energy conversion in the thermoelectric junction is usually characterized by the figure of merit 
(denoted as ZT), as defined in Eq. ([23|). ZT depends on the following physical factors: the Seebeck coefficient (S), 
the electric conductance (cr), the electron's thermal conductance (kez), and the phonon's thermal conductance (uph)- 
These physical factors can be evaluated using Eqs. ©, (IT5|) . p^. and (1551) . respectively. When ZT tends to infinity, 
the thermoelectric efficiency will reach the Carnot efficiency, the upper limit of energy conversion efficiency. To obtain 
a large ZT value, the thermoelectric nanojunction will require a large value of 5', a large value of ct, and a small 
value of the combined heat conductance (k = Kei + Uph)- Thermoelectric devices with a large value of a are usually 
accompanied by a large value of Kg;, due to the same proportionality with the transmission function. These values 
are highly correlated, making the enhancement of the thermoelectric figure of merit ZT a challenging task. 

The thermal energy carried by phonons flows from the hot into cold reservoir and is a negative effect to the 
thermoelectric refrigeration. We consider the rate of thermal energy carried by the phonon in the weak link model 
[as described in Eq. (jl9[) ]. where each electrode is assumed to be in thermodynamic equilibrium, joined by a weak 
mechanical link modeled by a harmonic spring of stiffness K. In the weak-tunneling limit, the weak link model can 
be interpreted as an application of the thermal Landauer formula [60| . Up to the leading order in temperatures, 
we expand Eq. ([T^ which gives a simple form, — X{T^ — T^) [i.e., Eq. ([551) specified as the simplified weak 
link model], where A w 2tt^ K'^ C'^ kg / {15h) . Simplification allows us to investigate the effect of the phonon's thermal 
current with a single parameter A, which is determined by the stiffness of the nano-structured object (K) and the 
slope of the surface phonon's spectral density (C) in the electrodes. The simplified weak link model is valid for 
T ^ Td, where To — 394 K is the Debye temperature for Al. We note that the simplified weak link model violates 
the Wiedemann- Franz law in low-temperature regime, while Eq. (1191) restores the Wiedemann- Franz law at T 3> Tjj. 
The range of temperatures discussed in this study lies within the valid regime of the simplified weak link model. 
Therefore, we consider ZT as a function of temperatures in the presence of the phonon's thermal current of which is 
represented by various values of A, as shown in Fig. [H^d). The validity of the A values for linear atomic wires will be 
justified later. 

B. Thermoelectric Cooling Effect in the Absence of Phonon's Thermal Current 

In the bulk system, the efficiency of a thermoelectric refrigerator is usually suppressed by the overwhelmingly 
irreversible Joule heating. Local heating in the nanojunction is analogous to Joule heating in the bulk system. Local 
heating arises from energy exchanged between the transport electrons and the vibration of the bridging nano-object. 
The heat generated in the nano-object via electron-vibration interaction is dissipated into electrodes, and is thus a 
negative effect to thermoelectric refrigeration. The amount of heat generated in the nanojunction depends on several 
factors, including the strength of electron- vibration interaction, background temperatures, and inelastic electron mean 
free path [22]. The nano-object bridging the junction is formed by few atoms, and thus the dispersion relation of 
phonon is characterized by the lack of Goldstone mode. When normal coordinates are considered, the complex 
vibrations of nano-object connecting to heavy electrodes can be cast into a set of independent simple harmonic 
oscillators described as normal modes. These features allow us to avoid the overwhelmingly irreversible heating effect 
by considering the 4-Al atoms junction operating under a bias much smaller than the onset bias Vonset (~ 20 mV for 
4-Al junction [6l|), where eVonset is the energy of the lowest normal mode and where the local heating is significantly 
suppressed. 

The suppression of local heating due to the small probability for electron-vibration interactions comes from two 
perspectives. First, the distance between two electrodes is much smaller than the electron-phonon mean free path, 
and thus the probability for electrons with sufficient energy to excite the vibration of nano-structure is small. Second, 
the atomic junction is operated under a small bias {Vb <C Vonset), and thus the majority of electrons have energies 
smaller than the lowest normal mode vibration of the nano-object connecting to the electrodes as a whole system, 
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FIG. 3: (Color online) (a) jJJ vs. Vb for various values of Tb. in the absence of the phonon's thermal current, where AT = 1 K; 
(b) JiJ vs. Tr for various values of AT in the absence of the phonon's thermal current, where Vb = 6 mV; (c) Jq vs. Vb in 



the presence of the phonon's thermal current for various values of A, where AT = 1 K and Tu 
the presence of the phonon's thermal current for various values of A, where AT = 1 K and Vb ■ 



: 300 K; and (d) vs. Tr in 
: 6 mV. 



apart from a small fraction of electrons that are thermally excited with energy larger that [ijx, + eVonset and due to 
the broadening effect. 

The atomic refrigerator possesses another prominent advantage, the suppression of photon radiation. Photon 
radiation is proportional to the surface size of the metal contact according to Stefan's law. The reduced surface size 
of the metal electrodes in typical nanojunctions significantly suppresses the energy transport carried by the photon 
radiation. Consequently, the photovoltaic effect, which is also against the thermoelectric cooling, can be neglected. 

The suppression of local heating and photon radiation due to size reduction remarkably facilitates the possible 
thermoelectric cooling effect in the 4-Al junctions beyond the expectation of the conventional bulk thermoelectric 
device theory. Our calculations indicate that as a thermoelectric nano-refrigerator, the 4-Al junction could work at 
temperatures below 100 K, which is a great improvement in comparison with the vacuum diode. In this regard, low 
operation temperature is an excellent characteristic of electronic nano-refrigerators. 

Let us now attempt to investigate the 4-Al junction as a thermoelectric cooling device from the first-principles ap- 
proaches. For the present, we shall confine our attention to the simplest case, the one which neglects the phonon's ther- 
mal current. The phonon's thermal current will be considered in the next subsection. We also put forward an analytical 
theory to explain the numerical results. The analytical theory described below is general to any atomic/molecular 
thermoelectric junction as a thermoelectric nano-refrigerator. 

When external bias is applied, electrons fiow from the right to the left electrodes. The fiow of electrons carries 
not only the charge current but also the energy current. The thermal current carried by electrons [denoted as 

J^i fJ-R,Tji)], defined as the rate at which thermal energy fiows from the right (into the left) electrode, can 

be calculated from first-principles using Eq. ([7]) with the help of the wavefunctions obtained self-consistently in the 
DFT calculations. Note that + IVb — J^i [for details see Eq. ([M)) ] because of energy conservation. It implies that 
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FIG. 4: (Color online) (a) COP (?;) vs. Vb in the absence of the phonon's thermal current for various values of AT, where 
Tr = 200 K; (b) rj vs. Vb in the absence of the phonon's thermal current for various values of Tr, where AT = 1 K; (c) r] 
vs. Vb in the presence of the phonon's thermal current for various values of A, where Tr = 200 K and AT = 1 K; and (d) rj 
vs. Vb for various values of Tr in the presence of phonon's thermal current, where AT = 1 K and A = 10~^^ J/(sec-K*). The 
enlarged graph (inset) shows that the zeros of COP (i.e., V^j^'^'''^) increases as Tr increases. 



the thermoelectric refrigeration requires electric power IVb to remove thermal energy from the cold (right)reservoir 
(with rate J^) and reject waste thermal energy to the hot (left) reservoir (with rate J^)- The nano-refrigerator works 
when > 0. As noted in Eq. ([7]), is a function of Vb, AT, and Tr. At a given Vb and AT, the solution of 
J^{Tr) = defines the threshold operation temperature, denoted as Top, below which the nano-refrigerator does 
not function. Figure [2^a) shows the Top of the 4 — Al nano-refrigerator as a function of AT and Vb- The operating 
temperature of the 4 — Al nano-refrigerator can be lower than 100 K, as shown in Fig.[2ja). 

To enrich the understanding of thermoelectric nano-refrigerators, we also propose analytical theories that provide 
guidelines for the design of nano-refrigerators. When AT <C Tr, the application of the Sommerfeld expansion on 
Eq. ([7]) can simplify the electron's thermal current as a polynomial of Vb, as described in Eq. (I?7|) . If the higher- 
order term Vg is neglected, the maximum of J^{Tji) that is greater than zero yields a criterion for the existence 
of thermoelectric cooling, which gives Eq. (UHl): -S > 7r/cB^2AT/(3Tc)/e, where Tc (= Tr) is the temperature of 
the cold (right) reservoir, S = — 7r^/c^T'(/i)Ti^/[3er(/i)] is the Seebeck coefficient, kB is the Boltzmann constant, and 
e is the electron charge. To have possible refrigeration effect, the Seebeck coefficient S, the temperature difference 
AT, and the temperature of the cold reservoir Tc need to satisfy Eq. (|28|) . The choice of an n-type thermoelectric 
junction with a large Seebeck coefficient is of crucial importance in the design of thermoelectric nano-refrigerators. 
Large Seebeck coefficients could be achieved thr oug h an appropriate choice of bridging nano-structured objects and 
further optimized by applying gate voltages [H, l48l | 

Figure [3ja) shows the rate of thermal energy extracted from the cold temperature reservoir (J^) as a function of 
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Vb for Tfj = 100, 200, and 300 K, where the temperature difference between the hot and cold reservoir AT is fixed 
at 1 K. We note that < at zero bias[see Fig. [3][a)], and thus a threshold voltage V^f for the battery is needed 
to trigger the refrigeration effect. The threshold bias is defined as the smallest positive solution of J^{Vb) = for a 
given Tr and AT. When Vb < V'j^', the rate of thermal energy removed from the cold reservoir is neg ative {J^ < 0) 
and thus the thermoelectric nano-refrigerator does not function. Since V^^' has a very small value, we neglect the 
terms higher than the first order of Vb in Eq. ([27| and obtain Vfil from J^{Vb) — 0, as described in Eq. ([29| . where 
we assumed AT ^ T^. Equation (|29p predicts that the threshold voltage Vf'jl — -~TT^kgAT/{3e'^S) is independent of 
Tr, as shown in Fig. Efa). 

Figure[3Ub) shows as a function of Tr for different values of AT — 0, 1, 2, 3, and 4 K, where the bias is fixed at 
Vb = Q mV. For a given Vb and AT, we define the critical operation temperature as the solution of J^{Tb) = 0. The 
thermoelectric nano-refrigerator is working when Tr > Top. The critical operation temperature can be calculated from 
Eq. (HID: Top « i[(a - 1)AT + y/{a^ - 1)AT^ + 13], where a = t{^il)/[t' iliL)eVB] and /? = -2Vb/S. The foregoing 
equation predicts that Top increases as AT increases, which agrees well with the numerical calculation presented in 
Fig. [Hb). If AT = 0, then Top reaches the minimum value (Top)^-^ « ^/TrV^/{^^2S) as described in Eq. The 
equation shows that {Top)^-^^ increases as Vb increases and {Top)^^^^ increases as T^ increases, respectively. At Vb = 6 
mV, the (2^op)min fo^ the i — Al thermoelectric refrigerator is 116 K, which can be further suppressed by decreasing 
the bias Vb ■ 

Turning to the coefficient of performance of the thermoelectric nano-refrigerator, the COP (denoted as 77*^') is defined 
as the ratio of the rate of heat removed from the cold reservoir to the electric power supplied by the battery, that is, 
77*^' — J^/IVb, as equivalent to Eq. ([5^ . Derived from Eq. (155)) . the maximum value of COP is given by Eq. ([Ml) . 
Figure |4l^a) shows the numeric calculations of 77*^' as a function of Vb for AT = 1, 2, 3, and 4 K, where Tr = 200 K. 
Figure lU^b) shows the numeric calculations of 77'^' as a function of Vb for Tr — 150, 200, 250, and 300 K, where 
AT = 1 K. The exchange of energy by heat currents between the hot and cold reservoirs is an irreversible process, 
and leads the decrease of the COP as AT increases, as shown in Fig. IH^a). This also agrees well with the prediction 
of Eq. dMl)- 

The maximum value of COP (denoted as t]^^^) occurs at bias Vj^^^^ ~ '^^trn where V^^ is the threshold bias for 
possible refrigeration given by Eq. (f^ : Ty^^x — 3e^5'^TR/(47r^fc^AT) — 1/2. As a direct consequence of the above 
equation, the maximum value of COP (?7max) increases as S'^Tr increases and decreases as AT increases. Remarkably, 
Vj2ax is independent of Tr as shown in Figs.lU^a) andlU^b). This result also agrees well with the theoretical prediction 
because Vj^ax ~ 2VJ^' and V^^ is independent of Tr as described in Eq. ([29]). We should note that the maximum 
value of COP can be greatly enhanced by a suitable nanojunction with a large Seebeck coefficient according to 
'7max ^ ^"^ predicted in Eq. (j34[) . Our study demonstrates that by taking advantage of the quantum features of 
reduced dimensions, a nano-refrigerator with a figure of merit ZT k, 0.02 has a COP comparable to a conventional 
TE refrigerator with a ZT « 0.2. 

C. Thermoelectric Cooling Effect Including the Phonon's Thermal Current 

The phonon's thermal current carries thermal energy from the hot to cold reservoir, which is an adverse effect 
to refrigeration. To assess the extent of this adverse effect, we consider the phonon's thermal current within the 
simplified weak link model. The simplified weak link model is suitable for describing the heat transport for two 
thermal reservoirs connected by a weak elastic link in a low temperature regime (T Tp) and is convenient to 
describe the phonon's thermal current by a single parameter A: — A(T^ — T^), where A w 2TT^K^C^k'^/{\bh). 
The parameter A can be determined by K and C, where K is the stiffness of the nano-structured object and C is 
the slope of the spectral density N{E) w C x E. From recent experiments (G^I, it is observed that the stiffness of 
the linear Pt atomic chains varies in a wide range of several orders of magnitudes. This experiment suggests that the 
stiffness K is likely to depend strongly on the detailed atomic structure of the full system, especially in the contact 
region. For example, the atom chain could be particularly stiff along the chain direction when the atomic chain forms 
a perfectly straight line; otherwise, the atomic (zigzag) chain could easily bend with much smaller stiffness. In this 
view, the magnitude of phonon's thermal current could possibly exist in a wide range of magnitudes in a nanojunction 
formed by a monatomic chain. An estimation from the experimental data of K and C in the Pt atomic junction 
yields the values of A range to be from to 2.05 x 10~^^ W/K'' in a typical atomic chain [gl, [6^. The simplified 
weak link model allows the construction of an analytical theory that offers a concise explanation for how the phonon's 
thermal current (represented by a single parameter A) affects thermoelectric cooling. We choose A, which may not be 
underestimated, around the order of 10~^^ W/K^ to compare with the value of Pt atomic junction. An alternative 
reason for choosing the A values around the order of 10^^^ W/K"* to present the phonon's thermal current is that 
thermoelectric refrigeration is sensitive with A around the order of 10^^^ W/K^ as shown below. 
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The COP in the presence of the phonon's thermal current J^^ is defined as 77 — jy^, where Jq = J^i + is 
the combined thermal current. In the presence of the phonon's thermal current, the functioning of the 4-Al atomic 
refrigerator requires higher operation temperatures, as shown in Fig. [2jb), where we plot the operation temperature 
Top as a function of Vb and AT with A = 10~^^ W/K^. Applying the Sommerfeld expansion, the combined thermal 
current Jq can be expressed as a polynomial of Vb when AT <IC Tr, as described in Eq. ([39|) . When the term 
Vg is neglected in Eq. ([39]). the criterion for the existence of electronic cooling is given by Eq. (|40l) . derived from 
the condition that the maximum value of Jq is positive. The criterion for possible thermoelectric refrigeration 
[Equation (|40p ] provides the guideline for devising experiments to test the effect of thermoelectric cooling at the 
atomic level. When A ^> 0, Eq. PO]) restores Eq. ([25]) in the absence of the phonon's thermal current. 

Figure [3ljc) shows the combined thermal current Jq as a function of Vb for A = 0, 3 x 10~^^ W/K^, and 6 x 
10"^^ W/K^, where we fix AT = 1 K and Tr = 300 K. Note that < at zero bias. Similar to the previous 
discussions in the absence of the phonon's thermal current, the thermoelectric nano-refrigerator works when Jq > 0. 

This leads to a threshold vohage [V^^^''^ = Vg + A [4 {~S) T^AT/a] , as given in Eq. (j4ip ]. which is the minimum 
external bias required to trigger the thermoelectric refrigeration. Equation (|41|) predicts that the threshold bias 
{Vg'^^^) increases as A increases, which agrees well with the numerical calculations as shown in Fig. ^c). Figure 
^d) shows the combined thermal current as a function of Tr for A = 0, 3 x 10"^'^ W/K", and 6 x 10"^'^ W/K^, 
where we fix Vb = 6 mV and AT = 1 K. Figure [3l[d) exhibits the combined thermal current Jq, which could become 
negative at high temperatures regime. The reason for this is that the phonon's thermal current brings heat from the 
hot to the cold reservoir. Thus, the combined thermal current becomes a negative value, which disables refrigeration 
capability at high temperatures. This imposes an upper limit for the working temperature, above which Jq turns to 
be negative and consequently the nano-refrigerator does not function. 

Derived from Eq. ([M)) . the upper hmit of COP (denoted as 'q'^tx'^ ) given by Eq. Vrntl'^ — 

r -1 1^1 

('ymax) + ^ (l6TRAT/(S'2cr)) . Equation (gS]) predicts that rjfJ^+P'' decreases as TrAT increases and <+f'' 

increases as S^a decreases. In Fig. [Ifc), we show the numeric calculations of the COP (jy'^'+P'*) as a function of Vb in 
the presence of the phonon's thermal current represented by A = 0, 1 x 10"^'', 3 x 10~^^, 6 x 10"^'' W/K^, where we 
fix AT = 1 K and Tr = 200 K. Figure HJc) shows that the upper hmit of COP {rj^tl'' ) decreases as the strength of 
phonon's thermal current increases because the phonon's thermal current is an adverse effect to refrigeration. Figure 
mjc) also shows that V^^i^~^^^ increases as A increases and agrees well with the prediction of Eq. (|4ip . Figure Sl^d) shows 
the numeric calculations of the COP (ry«='+P'') as a function of Vb for Tr — 200, 250, 300, and 350 K, where we fix 
AT = 1 K and A = 1 x 10^^'' W/K'^. As shown in Fig. |4l^d), the upper limit of COP (?7max) increases as A increases as 
predicted by Eq. ([1^ : meanwhile, the inset of Fig. HJd) shows that the Vg~^^'^ increases as Tr increases, as predicted 
by Eq. (gH) when AT < Tr, 

The phonon's thermal current, which is relevant to the mechanical coupling between the nano-structured object 
and the electrodes, is an adverse effect to refrigeration. To minimize the adverse effects, we suggest creating a weak 
mechanical link between the nano-structured object and the electrodes while still allowing electrons to tunnel. 



IV. CONCLUSIONS 



In theory, we propose an innovative idea for finding a possible thermoelectric cooling device at the atomic scale. 
In-depth research from the first-principles approaches is performed to investigate the thermoelectric cooling effect of 
the atomic-scale refrigerator. This research offers a key to understanding the energy-conversion mechanism at the 
atomic scale. Our studies show that the choice of thermoelectric junction with a large Seebeck coefhcient is of crucial 
importance for the design of thermoelectric nano-devices. The bridging Al atoms create potential wells that lead 
to resonant tunneling via a sigma channel near the chemical potentials with the Px-Py orbital characters, and thus 
significantly enhance the Seebeck coefficient of the 4-Al atomic junction. 

When the nano-refrigerator is operated under a bias smaller than the onset voltage (Vb ^ Vonset), the local heating 
is significantly suppressed. Our research shows that the suppression of local heating, photon radiation and tunneling 
barrier by the bridging atoms in nanojunction is favorable for the operation of thermoelectric nano-refrigerators. 
The meritorious factors mentioned above remarkably facilitate thermoelectric cooling beyond the expectation of the 
conventional thermoelectric device theory in the bulk system. 

Analytical theories have been established to offer a more elaborate explanation of the working conditions, threshold 
bias, and properties of coefficient of performance. We observe a threshold bias Vth below which no refrigeration is 
possible. The COP of nano-refrigerators is optimized at Vb ~ "^^tii when AT ^ Tr, A major adversary factor to 
thermoelectric refrigeration is the phonon's thermal current, which is likely to depend strongly on the detailed atomic 
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structure of the full system, especially in the contact region. To minimize the adverse effect of the phonon's thermal 
current, we suggest creating a failure in the mechanic link between the nano-structurcd object and the electrodes, 
while still allowing the electrons to tunnel. 

This research has theoretically revealed that an atomic-scale thermoelectric refrigerator is promising. The nano- 
refrigerators potentially have COP comparable to conventional TE refrigerators with ZT, which is an order of mag- 
nitude larger. Moreover, widely diversified atomic-sized junctions may be achieved by manipulating the species of 
nano-structured objects and the contact region. Such manipulations may lead to a significant change in the DOSs 
and the stiffness of nanojunction, consequently varying the Seebeck coefficients and the phonon's thermal current. 
Atomic-level control of the contact region is expected to open new opportunities and challenges in developing new 
forms of thermoelectric energy conversion devices, which need to be extended by the utilization of unprecedented 
experiments. We hope that our forward-looking theoretical proposal could stimulate innovated experiment towards 
further exploration of atomic / molecular system which may lead to practical applications of nanoscale thermoelectric 
cooling devices in the future. 
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